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ABSTRACT 

We have studied the sensitivity of s-process nucleosynthesis in massive stars to ±2(7 variations in 
the rates of the triple-a and ^^C(a,7)^^0 reactions. We simulated the evolution of massive stars 
from H-burning through Fe-core collapse, followed by a supernova explosion. We found that: the 
production factors of s-process nuclides between ^^Fe and ^^Zr change strongly with changes in the 
He burning reaction rates; using the Lodders (2003) solar abundances rather than those of Anders 
& Grevesse (1989) reduces s-process nucleosynthesis; later burning phases beyond core He burning 
and shell C burning have a significant effect on post-explosive production factors. We also discuss the 
implications of the uncertainties in the helium burning rates for evidence of a new primary neutron 
capture process (LEPP) in massive stars. 

Subject headings: nuclear reactions, nucleosynthesis, abundances, Supernovae: General 



1. INTRODUCTION 

Half of the elements between Fe and Bi are produced 
by the slow (s) neutron capture process and most of 
the remainder by the rapid (r) neutron capture pro- 
cess. About 35 additional neutron deficient stable nu- 
clides above ^^Fe, the p-process nuclei, are produced in 
explosive processes (Prantzos et al. 1990). 

Two components, the main and the weak s-processes, 
are required to explain the isotopic distributions of s- 
process nuclides. The main s-process occurs in low-mass 
(< 4M0) asymptotic giant branch (AGB) stars and con- 
tributes mainly to the production of heavier elements, 
with a smaller contribution to ^4 < 90. The weak s- 
process occurs during the late evolutionary stages of mas- 
sive stars (^10 M0) and produces nuclides up to A :^ 90. 

Recently, Travaglio et al. (2004) summed the contribu- 
tions of the r-process, the main and weak s-process, and 
the p-process to the abundances of Sr, Y, and Zr. The 
summed contributions were smaller than the observed 
solar abundances by 8%, 18%, and 18%, respectively. 
They concluded that an additional light element pri- 
mary s-process contribution from massive stars (LEPP) 
is needed to explain this difference; the nature and site 
of the LEPP are unknown. This LEPP has also been in- 
voked by Montes et al. (2007) to explain the abundances 
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of a larger group of light r-process elements. Since the 
LEPP effects are relatively small, and since some of the 
LEPP elements are produced with relatively large abun- 
dance in the weak s-process, it is important to establish 
whether the uncertainties in the weak s-process are suf- 
ficiently small that the claim of LEPP contributions is 
robust. 

Nuclide production in the weak s-process also depends 
on the rate of the neutron source ^^Ne(a, n)^^Mg and 
on the capture cross sections for the neutron poisons 
(medium- weight isotopes up to Fe, including ^^C, ^^O, 
2^Ne, ^^Ne, ^^Na, ^"^Mg, and ^^Mg). Neither the source 
strength nor the neutron capture cross sections for the 
poisons are known with sufficient accuracy (Heil et al. 
2008). We shall not deal with these issues in this paper, 
but rather with the more indirect effects of uncertain- 
ties in the rates {Rsa and Ra,i2) of the triple alpha and 
^^C(a,7)^^0 reactions, and in the initial stellar compo- 
sition. For example, we have shown in a previous pa- 
per (Tur et al. 2007; see also Weaver & Woosley 1993; 
Woosley et al. 2003; Woosley & Heger 2007), that the 
amount of the above neutron poisons present during the 
weak s-process in massive stars depends sensitively on 
these rates and the initial stellar composition. 

Most earlier studies of the weak s-process focused on 
production toward the end of core He burning by neu- 
trons from the ^^Ne(a, n)^^Mg reaction (Couch et al. 
1974; Lamb et al. 1977; Arnett et al. 1985; Busso & 
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Gallino 1985; Langer et al. 1989; Prantzos et al. 1990; 
The et al. 2000; Raiteri et al. 1991a; Baraffe et al. 1992). 
Later papers considered also a second exposure at higher 
temperatures and neutron densities peaking during shell 
C burning (Raiteri et al. 1991b; Raiteri et al. 1992; Rai- 
teri et al. 1993; The et al. 2007). Explosive processing 
in the supernova explosion was not considered. Recent 
calculations of s-process yields have been extended to 
consider the entire evolutionary history of the star, in- 
cluding explosive burning (Hoffman et al. 2001; Rauscher 
et al. 2002; Limongi & Chieffi 2003). 

In this paper we simulate the evolution of massive stars 
from H-burning through Fe-core collapse, followed by a 
supernova explosion. Our principal purposes are (1) to 
establish the magnitude of weak s-process production in 
a self-consistent model; (2) to study the uncertainties in 
weak s-process nucleosynthesis arising from uncertainties 
in Rsa and i^a,i2; (3) to study the effects of different stel- 
lar abundances, specifically those of Lodders (2003) and 
Anders & Grevesse (1989), hereafter L03 and AG89; (4) 
to delineate the stages of stellar evolution during which 
weak s-process elements are produced; and (5) to assess 
the bearing of these uncertainties on the robustness of 
the LEPP process. 

In Section 2, we describe the stellar model and the 
input physics relevant to the treatment of the weak s- 
process. Section 3 gives our results for the dependence of 
the post-explosive weak-s process production factors on 
changes in the rates of the helium burning reactions and 
on different initial stellar abundances. In Section 4, we 
show the contribution to the weak s-process abundances 
of the various stellar burning stages prior to the super- 
nova explosion. In Section 5, we investigate the range of 
weak s-process production of Sr, Y, and Zr allowed by 
the uncertainties in the helium burning reactions. 

2. STELLAR MODELS AND INPUT PHYSICS 

The one-dimensional hydrodynamics code KEPLER 
(Weaver et al. 1978; Woosley & Weaver 1995; Rauscher et 
al. 2002; Woosley et al. 2002) is used to evolve stars from 
central H burning up to Fe core-collapse. The supernova 
explosion that follows is simulated using a spherical pis- 
ton placed at the base of the O burning shell, which 
first moves inward and then outward at a constant ac- 
celeration which has been adjusted to result in a total 
kinetic energy of the ejecta of 1.2 x 10^^ ergs one year 
after the explosion, and coming to rest at a radius of 
1,000km (Rauscher et al. 2002; Woosley & Heger 2007). 
After estimating the fallback from our hydrodynamic su- 
pernova simulations, the final nucleosynthesis yields are 
determined by employing the same parameterization of 
mixing as was used by Woosley & Heger (2007). 

We performed calculations for stars of 15, 20, and 
25 M0 and for two different initial abundances, L03 and 
AG89. Rate sets were: (1) Rsa was kept constant (at 
its value from Caughlan & Fowler 1988), and Ra,i2 was 
varied; (2) both rates were varied by the same factor, 
so their ratio remained constant; and (3) Ra,i2 was held 
constant at 1.2 times the rate recommended by Buch- 
mann (1996) and Rsa was varied. Both reaction rates 
were varied within a range of ±2cr of their experimental 
uncertainties. 

For more details see Rauscher et al. (2002); Tur et al. 
(2007); Woosley & Heger (2007) who give a complete 



description of the improvements to the stellar physics 
and reaction rates since Woosley & Weaver (1995). In 
particular, the "rath" rates, Rauscher & Thielemann 
(2000), have been adopted for the competing reactions 
22Ne(a,n)25Mg and ^^Ne{a,j)^^Mg. 

We adopt a two-character notation to label our plots, 
e.g., LA, similar to the notation adopted by Tur et al. 
(2007). The first character can be an L (to denote the 
L03 initial abundances) or an A (for the AG89 initial 
abundances). The second character denotes the study: 

A, when Rsa was kept constant, and Ra,i2 was varied, 

B, when both rates were varied by the same factor, so 
their ratio remained constant, and C, when Ra,i2 was 
held constant, and Rsa was varied. 

3. POST-EXPLOSIVE PRODUCTION FACTORS 

In Figure 1, we show the evolution of the production 
factors as a function of Ra,i2 for isotopes between ^^Fe 
and ^^Zr. To facilitate reading. Figure 2 shows the same 
for the 6 s-only isotopes. We define the production factor 
for each stable isotope as the ratio of the mass fraction 
of the isotope in the supernova ejecta to its initial solar 
mass fraction. Rsa was kept constant at the value from 
Caughlan & Fowler (1988). The color bar range is from 
1 to 3 X 10^. The range of production factors for a 25 Mq 
star are shown in Table 1. For most assumed reaction 
rates, production factors (PF) obtained using the L03 
abundances (Figures 16, Id^ and 1/), are significantly 
smaller than those for the AG89 abundances (Figures 1 a, 
Ic, and le). Presumably the lower CNO content of the 
L03 abundance set leads to a smaller amount of ^^Ne, 
lower neutron abundance, and a less efficient s-process 
(Woosley & Heger 2007). 

We have shown in Tur et al. (2007) that for the range of 
reaction rates we calculate, the PF for ^^O is 10.0 ± 0.5 
for the AG89 abundances, and 15.3 =b 0.5 for the L03 
abundances; the uncertainty is the standard deviation 
of the results for the different rates. We find that the 
PFs we obtain for the s-process elements are often larger 
than this: see Figure 1 and Table 1. This is not surpris- 
ing: very low metallicity stars produce oxygen, but very 
small abundances of the secondary weak s-process nuclei. 
We next use a simple model to estimate what PF in a 
solar-abundance star would be required to reproduce the 
observed abundances. 

Assume that a nuclide X is to be produced entirely 
by secondary processes in massive stars. It follows 
that when integrated over metallicities from zero to the 
present, its PF must be the same as that for ^^O. Solar 
metallicity stars must, therefore, have a larger PF for 
X than for ^^O. In a simple closed-box model, assuming 
the production of X is proportional to metallicity, the 
required ratio of secondary to primary production in a 
solar abundance star is two. We take this as a work- 
ing hypothesis, bearing in mind that it is only a rough 
approximation. Then, if the PF for a weak-s nuclide 
is 2xPF(^^0), massive stars produce the observed solar 
abundance, or if it is equal to PF(^^O), half the observed 
abundance. The fraction / of the observed solar abun- 
dance for nuclide X produced in our model is given by 

/(X) = 0.5 X PF(X)/PF(^^0) (1) 

where the 0.5 arises from the model described above; as 
noted above PF(^^O) is approximately 10 (15.3) for the 
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Fig. 1. — Post-explosive production factors as a function of Ra,i2- 
eries, 20 M©. (d) The LA series, 20 M©. (e) The AA series, 25 M© 
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(a) The AA series, ISM©, (b) The LA series, ISM©, 
(f) The LA series, 25 M©. 



(c) The AA 



AG89 (LOS) abundances. Equation (1) and Figure 1 
can be used to estimate whether a particular isotope is 
overproduced for a given star or rate choice. 

The production factors increase strongly as the stellar 
mass increases, reflecting a more efficient s-process as 
the temperature increases. In particular, most isotopes 
are weakly produced, PF<20 (30) for AG89 (LOS) in the 
15 M© star, with the exception of those between ^^Fe and 
^^Cu, and of some isotopes when Roc,i2 is given its lowest 



or highest value (Figures la and Ih). The highest PF are 
observed for the 25 M© star (Figures le and 1/), whereas 
the 20 M0 star shows a contribution between those two 
extremes (Figures Ic and Id). Ignoring the contribution 
of the 15 Mq and 20 stars as was done in Raiteri 
et al. (1991a) and Raiteri et al. (1991b) does not seem 
justified; their contributions are not always negligible. 

We find that ^^Ni is overproduced for many rate 
choices, especially for the 25 star and the AG89 
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^^C(a,7)^^0 rate multiplier ^^C(a,7)^^0 rate multiplier 

Fig. 2. — Post-explosive production factors for s-only isotopes as a function of Ra,i2- (sl) The AA series, 15 M©. (b) The LA series, 
15 M0. (c) The AA series, 20 M©. (d) The LA series, 20 M©. (e) The AA series, 25 M©. (f) The LA series, 25 M©. 



abundances. This problem was already documented by 
Rauscher et al. (2002) and seems to arise from the de- 
creased value of the ^^Ni(n, 7)^^Ni reaction rate used for 
our models compared to its previous compilation. The 
factor of three change arises from different extrapolations 
of a cross section measured only at thermal neutron en- 
ergies (see also Nassar et al. 2005, however). 

Table 1 shows, for a 25 star of L03 initial abun- 
dance, the value of the maximum and minimum produc- 



tion factors for nuclides between ^^Fe and ^^Zr when ei- 
ther Ra,i2 or Rsa is varied within its ±2a uncertainties. 
Large differences in PFs are found when either Rsa or 
Ra,i2 are varied. Most commonly, the largest PFs are 
obtained for the lower values of Ra, 12 as one sees in Fig- 
ure 1. Figure 3 shows some of the changes resulting from 
changes in Rsa- Lower values of Rsa lead to lower ^^C 
abundance at the end of central helium burning. This, in 
turn, is compensated by a slightly higher central temper- 
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0.00") for 25 M0 stars with L03 initial abundance as a 
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Fig. 5. — Post-explosive production factors for s-only isotopes for a 25 Mq star as a function of Rsa- (a) The LC series, (b) The LB 
series. The "reaction rate multiplier" is the factor that multiplies both the standard rates. 
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Fig. 6. — Energy history of the 25 star (central values of the 
reaction rates, L03 abundances) as a function of time until core 
collapse. The ordinate is the included mass from the stellar cen- 
ter. The green cross hatched areas are fully convective, and the red 
cross hatched areas are semiconvective. The blue and pink shading 
indicate net energy generation, with blue positive and pink nega- 
tive. For more details, see Woosley et al. (2002). The vertical lines 
show the times when dump files are generated by the KEPLER 
code at key stages of the evolutionary process. 



ature to produce the same energy release rate to main- 
tain the star's luminosity. Due to the high temperature 
sensitivity of the ^^Ne(a,n) reaction, this increased tem- 
perature means more of the ^^Ne has been burnt by cen- 
tral helium depletion, resulting in a stronger s-process. 
Figure 3 also shows the decrease of ^^Fe for lower Rs^ 
that correlates very well with the smaller amount of ^^Ne 
left, i.e., more of it being burnt, providing a measure of 
neutron exposure. 

For many isotopes, the maximum PF is greater than 
twice the minimum. For example, for ±2(7 uncertainties 
in Ra,i2 the PF for ^^Ge ranges from 1.5 to 6.5; 4.2 to 
20.9; 12.1 to 38.0 for 15 M©; 20 M®; 25 Mq stars, resp. 
The s-only isotopes show a similarly strong sensitivity. 
For ±2(7 uncertainties in i^a,i2, the PF for ^^Ge, an s- 
only isotope, decreases from 6.6 to 2.2; 18.8 to 9.5; 32.5 
to 10.9 for 15 Mq; 20 M©; 25 Mq stars. For an average 
over the 3 stars, the production factor for ^^Kr, another 
s-only isotope, ranges from 2.5 at the lowest value of 
Ra,i2 to 20.5 at its next to highest value. 

In Figure 4, we illustrate the variations in the produc- 
tion factors for a 25 star with the L03 initial abun- 
dances when Rsa is varied and Ra,i2 is held constant 
(Figure 4a), and when both reaction rates are varied 
by the same amount (Figure 46). To facilitate reading, 
we show the same for the 6 s-only isotopes in Figure 5. 
Large variations in the production factors are observed 
in both cases, demonstrating the importance of both he- 
lium burning reactions. 

The results shown in Figure 1 and Figure 4 are con- 
sistent with the findings of Tur et al. (2007) that both 
reactions sensitively infiuence the production factors of 
the medium- weight isotopes, among which are the neu- 
tron poisons mentioned earlier. When, however, we ex- 
amine the production of these poisons and of the weak-s 
isotopes, as a function of the rates, we find only weak 
indications of any correlation. 



In Table 2, we compare our results for the s-only el- 
ements, using the central values of the reaction rates, 
with results of earlier calculations. The temperature de- 
pendence of the ^^Se beta decay rate is not implemented 
in KEPLER; this may account for the lower value of the 
production factor for ^^Kr compared to the results of The 
et al. 2007. Our values of weak-s production percentage 
are significantly larger than those of Raiteri et al. (1993) 
and The et al. (2007) presumably because of contribu- 
tions from burning beyond O depletion. The present re- 
sults have the correct trend, being smaller when the main 
s-process is large. Overall, however, the results for these 
central values of the rates are surprisingly large. But 
they are also very sensitive to the reaction rates. For ex- 
ample, for the L03 abundances, a 15 % smaller value of 
the ^'^C{a,j)^^0 rate gives weak-s contributions at least 
a factor of 2 smaller for these nuclei. On the other hand, 
for the AG89 rates, the weak-s contributions would be 
still larger. This is yet another example of the sensitivity 
of supernova nucleosynthesis to the rates of the helium- 
burning reactions, and further evidence of the need for 
more accurate rates. 

4. EVOLUTION OF PRODUCTION FACTORS ALONG THE 
STELLAR BURNING HISTORY 

In this section, we compare the contribution of various 
stellar burning phases to the s-process abundances for 
the 25 M0 star, L03 initial abundances, and the central 
values of the helium burning reaction rates: Rsa from 
Caughlan & Fowler (1988), and Ra,i2 equal to 1.2 times 
the rate recommended by Buchmann (1996). Dump files 
are generated at specified times during the life of the 
star as shown in Figure 6. From these we extract the 
PFs for isotopes between ^^Fe and ^^Zr and plot them 
in Figure 7 versus their distance from the stellar center 
in solar masses. The PF is defined as the mass fraction 
of an isotope in a given stellar zone (or mass location) 
divided by its solar mass fraction. All the color bars of 
Figures 7a-h have the same scale (from 1 to 4 x 10^) to fa- 
cilitate comparison among them. Figures 7a-g show the 
evolution of production factors at various burning phases 
of the star, up to the pre-supernova stage. The black ver- 
tical line on the figures shows the location of the initial 
mass cut (i.e., the position of the piston at the base of 
the O burning shell right before the explosion). Only the 
stellar mass shells above that line have a chance of being 
ejected during the explosion. The gray shaded areas on 
figures 7/ and 7g show those regions of the star where nu- 
cleosynthesis calculations are no longer being performed 
because a nuclear statistical equilibrium network is em- 
ployed and these layers are known to become part of the 
iron core and will not mix with the layers above before 
core collapse. They also lie inside the piston, i.e., below 
the pre-supernova mass cut. 

Below we give a detailed explanation for Figure 7: 

• Fig. 7a: at He ignition. Since no s-processing has 
yet occurred at this stage, the production factors 
are essentially zero. 

• Fig. 7b: at He depletion (when the central He 
mass fraction has reached 1 %). Central helium s- 
processing is in progress during this late He burning 
phase, as the temperatures become high enough to 
ignite ^^Ne(a, n)^^Mg. 
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Fig. 7. — Production factors for s-process isotopes for a 25 star (L03 initial abundances) as a function of their mass location within 
the star, (a) At He ignition, (b) At He depletion (1% He left), (c) After central He depletion but before C ignition, (d) At central 
C depletion, (e) At central O depletion, (f) At central Si depletion, (g) At onset of core collapse, (h) 100 s after core collapse (after 
explosive nucleosynthesis). 
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Fig. 8. — Production factor versus the mass coordinate (1.2 — S.SM©) for a 25 star viewed at various times during the star's hfe 
for (a) ^'-'Ge (b) ^^As (c) ^'-'Kr (d) ^^Rb. The figures are plotted for the central values of the helium burning reaction rates. The thick 
vertical line shows the location of the initial mass cut. 
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• Fig. 7c: just before central C ignition (when the 
central temperature reaches a value of 5 x 10^ K). 
Central helium burning has finished; this figure 
shows the final abundance and distribution of iso- 
topes due to central helium burning. About 1/3 of 
this s-process yield is located inside the later mass 
cut (compact remnant) of the supernova. This has 
been largely neglected in most of the previous pa- 
pers on the subject as it is hard to assess without 
simulation of the supernova. 

• Fig. 7d: central C depletion (when the central 
temperature has reached 1.2 x 10^ K). Central C 
burning contributes mostly to those regions of the 
star below the mass cut. We also see the result 
of onset of the first carbon burning shell and of 
He shell burning (between mass coordinate of GM© 
and IMq) which appears to be significant for most 
isotopes. 

• Fig. 7e: at central O depletion (when the central O 
mass fraction has dropped below 5%). Significant 
s-processing occurs during the C shell burning as 
seen in the enhancements to the PFs in the mass 
region above the mass cut. There is also photodis- 
sociation during central O burning in the region 
below the mass cut. 

• Fig. 7/: at Si depletion fwhen the central Si mass 
fraction drops below 10~ ). Additional photodisso- 
ciation occurs during core Si burning. Very strong 
overproduction occurs for some isotopes close to 
the mass cut. 

• Fig. 7g: at the pre-supernova stage (when the 
contraction speed in the iron core reaches 1000 
kms""*^). 

• Fig. 7h: 100 s after the passage of the shock wave. 
This panel shows the PFs above the mass cut, when 
most nuclear reactions cease due to the very low 
temperatures. The isotopes closest to the mass 
cut have essentially been destroyed with the excep- 
tion of some local overproduction for isotopes be- 
low ^^Ge. These local enhancements can still lead 
to over-all increased supernova yields after mixing 
and fallback. 

In Figures 8a-(i, we show the production factors of sin- 
gle isotopes (^^Ge, ^^As, ^^Kr, and ^^Rb) for a 25 M© 
star versus the mass coordinate (1.2 — 3.5 M©) at var- 
ious times during the star's life, using the central values 
of the helium burning reaction rates. Figure 8 c highlights 
the importance of the late evolutionary stages: most of 
the ^^Kr made during the stages up to pre-supernova 
is destroyed and rebuilt during the supernova explosion. 
For '^^Ge, '^^As, and ^^Rb, the C sheh burning and the 
pre-supernova stages both have significant contributions 
for masses beyond about 2.5 M©. The temperature de- 
pendence of the ^^Se beta decay rate is not implemented 
in KEPLER; the low values of the production factor for 
^^Kr in the earlier evolutionary stages presumably re- 
flect this fact. This comment also applies to Figure 9. 
The lower values of the ^^Ne(a, n)^^Mg rates used here 
may also contribute to smaller yields compared to some 
previous studies. 
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He-dep C-ign C-dep 0-dep Si-dep pre-SN Final 
STAGES 

Fig. 9. — Production factors of s-only nuclei lying outside the 
mass cut at various evolutionary stages for a 25M0 star using the 
central values of the reaction rates and the L03 abundances. 



To assess the importance of various stages we have de- 
termined the PFs for isotopes produced outside the mass 
cut at the times shown in Figure 6. It is clear from Fig- 
ure 9 that the later evolutionary stages contribute sig- 
nificantly, especially for ^^Kr. In the cases of ^^'^^Sr, 
however, little changes after oxygen depletion. 

5. SR, Y, ZR, AND THE LEPP PROCESS 

As noted in the introduction, the LEPP process was in- 
troduced by Travaglio et al. (2004) to explain the excess 
of the observed abundances of some elements, especially 
Sr, Y, and Zr, above that produced by the r-process, the 
main and weak s-processes, and the p-process. We exam- 
ine here the uncertainties in the weak s-process produc- 
tion of these elements. The differences shown in Table 1 
and Fig. 1 for individual isotopes appear much larger 
than the 8 % — 18 % differences that led to the proposal 
of the LEPP process. These differences, however, are 
important only if, (1), the weak s-process makes a suf- 
ficient fraction of an isotope, and if, (2), the differences 
persist for the elemental differences, the observed quanti- 
ties, not just the differences in a single isotope. We have 
calculated the elemental production factors for the 25 Mq 
star and the integrated result for the 15 Mq, 20 Mq and 
25 Mq stars, using the nominal reaction rates for R^^a 
and R(x,i2 and the LOS abundances. The results were 
averaged over an IMF with a slope of 7 = —2.6 (Scalo 
1986) and divided by their solar mass fraction, as de- 
scribed in Tur et al. (2007). We summarize the relevant 
quantities in Table 3. For the 25 Mq star taken alone, 
the uncertainties are significant, 10 % — 16 % compared 
to the 8 % — 18 % differences that were to be explained 
by the LEPP process. When one takes the averages over 
the 15 Mq, 20 Mq, and 25 Mq stars, however, uncertain- 
ties are smaller, 3% — 7%, but those for Sr are similar 
to the 8 % LEPP effect. All of these uncertainties are 
upper limits, in the context of the present calculations, 
since they are the extreme values for the entire set of re- 
action rates. Differences in PFs for these elements tend 
to cancel when one adds the production for the three 
stars, but it is not known whether such cancellations will 
occur for other elements. And one must still consider 
the uncertainties arising from uncertainties in the rates 
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of the neutron producing cross sections and the capture 
cross sections for various poisons as described in the in- 
troduction. Taken together, these considerations leave 
the LEPP process only moderately robust. 

6. CONCLUSIONS 

We followed the entire nucleosynthesis throughout the 
life of massive stars, from H-burning through Fe-core col- 
lapse, followed by a supernova explosion. We observe a 
strong sensitivity of the PFs for weak-s process isotopes 
to ±2<j variations of the rates of the triple-a and of the 
^^C(a,7)^^0 reactions. This can be explained by the 
variations in the PFs of medium-weight neutron poisons 
found by Tur et al. (2007), the changes in the amount of 
carbon left at the end of central helium burning, and the 
amount of ^^Ne burnt in central helium burning. 

In most cases, our simulations yield lower PFs for the 
abundances of LOS than for the abundances of AG89; 
the lower CNO content of the LOS abundance set is re- 
sponsible for the reduced efficiency of the s-processing. 
This tendency is not always followed, however, as we see 
for the s-only nuclei with the central rates, in Table 2. 
The production of the weak-s nuclei is highly sensitive 
to the rates of the helium burning reactions. In some 
cases (see the discussion of the results of Table 2), we 
find that a 15 % change in these rates may change the 
nucleosynthesis rates by more than a factor of two. 

We find that one must follow the entire evolution of the 
star to evaluate accurately the contribution of massive 
stars to s-process abundances. Most earlier studies took 
into account only the contribution of core He burning 
and shell C burning. We show that significant produc- 
tion takes place in later burning phases as detailed for 
the s-only isotopes in Figure 9. The burning phases be- 
yond shell C burning mostly destroy isotopes in the core 
(especially near the mass cut), but overall lead to an in- 
creased production of some isotopes. The passage of the 
shock wave further modifies those PFs, usually reducing 
them slightly. 

We have examined the uncertainties in the weak-s pro- 
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cess production for Sr, Y, Zr owing to uncertainties in the 
heUum-burning reaction rates, and find that they are in 
the 3 % — 7 % range. This is smaller than the observed 
deficiencies of 8 % — 18 % in the known production mech- 
anisms that led to the introduction of the LEPP process. 
On the other hand, we surely underestimate the total 
uncertainties. Uncertainties in the reaction rates of the 
neutron producing and capture reactions are significant. 
In addition, uncertainties in the stellar models can have 
significant effects (e.g., the treatment of hydrodynamics, 
convection, overshoot, etc.; see Costa et al. 2007 for a 
study of the effects of overshooting). We refer the reader 
to Tur et al. (2007) for a more complete discussion of 
those approximations and physics uncertainties. Com- 
bining all these uncertainties may render the evidence 
for the LEPP process less convincing. 

We have shown that many aspects of nucleosynthesis in 
supernovae, the production of the medium weight nuclei, 
as discussed in Tur et al. (2007), and that of the weak- 
s-process nuclei described here, are highly sensitive to 
variations of the helium burning reaction rates, within 
their experimental uncertainties. This further empha- 
sizes the need for better values of the helium burning 
reaction rates. 
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Table 1. Minimum and maximum values of the s-process production 
factors for a 25 m0 star with los initial abundances 



Isotope 


Within 2cr 


errors in Ra,i2 


Within 2a 


errors in Rsa 


Maximum 


Minimum 


Maximum 


Minimum 


58 pg 


18.8 


13.7 


19.1 


13.1 


^9 Co 


15.6 


9.3 


12.0 


9.2 




16.2 


8.7 


11.9 


6.0 


6lNi 


36.1 


20.7 


29.6 


20.7 




52.5 


22.8 


31.3 


22.8 




43.8 


21.0 


37.7 


21.0 


^^Cu 


35.0 


14.1 


33.2 


14.3 




34.5 


18.5 


26.5 


19.7 


64 Zn 


4.4 


2.3 


3.7 


2.2 


66 Zn 


19.9 


10.8 


20.0 


10.5 


67Zn 


26.3 


11.4 


22.9 


11.4 


68 Zn 


22.8 


12.2 


20.0 


12.2 


69Ga 


36.0 


17.3 


32.6 


21.8 




43.8 


19.2 


40.7 


19.2 


'^^Ge 


32.5 


10.9 


36.3 


13.8 


'^^Ge 


23.2 


12.4 


28.7 


11.9 


'^^Ge 


38.0 


12.1 


30.2 


11.7 


'^^Ge 


25.1 


8.0 


16.8 


8.0 




30.4 


13.4 


30.2 


12.9 


'^6Se 


43.5 


17.5 


54.3 


13.9 




35.8 


12.2 


30.3 


12.2 




15.6 


7.6 


20.0 


6.4 


80Se 


13.3 


5.3 


10.4 


5.3 


79 b J. 


22.7 


7.2 


17.8 


9.3 




13.8 


8.3 


16.7 


6.7 


80 Kr 


42.6 


4.7 


47.6 


4.1 


82 Kr 


26.2 


7.3 


29.9 


5.5 


83 Kr 


12.7 


5.3 


12.6 


5.3 


84 Kr 


6.0 


4.1 


6.6 


3.8 


86 Kr 


48.4 


3.1 


13.1 


3.1 


85 Rb 


19.5 


6.8 


16.9 


6.8 


87Rb 


20.5 


7.4 


18.8 


6.4 


86 Sr 


20.3 


3.7 


23.9 


2.9 


87Sr 


15.3 


2.3 


13.1 


2.9 


88 Sr 


6.4 


3.4 


5.8 


3.3 


89y 


4.7 


2.6 


5.8 


2.8 


90Zr 


6.8 


2.1 


8.0 


1.9 


9iZr 


3.1 


1.6 


3.9 


2.1 


92Zr 


3.1 


1.5 


3.2 


1.5 


94Zr 


2.4 


1.2 


2.0 


1.2 


96Zr 


5.3 


1.0 


3.3 


1.0 



Table 2. Percentage of s-only nuclei produced by the weak s-process 
for a 25 m0 star 





PF(X) AG 


PF(X) The^ 


PF(X) Lod 


PF(AG)/PF(Lod) 


PF(X)/PF(i6o)AG 


PF(X)/PF(i6o)Lod 


%Weak-s AG^ 


%Weak-s Lod^ 


%Weak-s Raiteri^ 


70 Ge 


24.9 


32.9 


32.5 


0.8 


2.49 


2.12 


124.7 


106 


64 


76Se 


38.7 


29.2 


39.4 


1.0 


3.87 


2.58 


193.6 


129 


63 


80 Kr 


18.5 


47.7 


28.8 


0.6 


1.85 


1.88 


92.7 


94 


86 


82 Kr 


20.9 


28.4 


20.0 


1.0 


2.09 


1.31 


104.5 


65 


53 


86 sr 


7.5 


19.9 


18.1 


0.4 


0.75 


1.18 


37.4 


59 


24 


87Sr 


9.1 


15.8 


10.1 


0.9 


0.91 


0.66 


45.7 


33 


16 



^From The et al. (2007); results labelled 25K in Table 7. 

^Calculated using Eq. 1. 

^From Raiteri et al. (1993); Table 5. 
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Table 3. Fraction of Sr, Y, Zr produced in the weak-s process LOS 

ABUNDANCES 



PF-max PF-min f-max^ f-min^ Af^ 



Sr 8.17 
Y 5.76 
Zr 5.24 



25 M0 star: 

3.30 0.27 
2.64 0.19 
2.09 0.17 



0.11 0.16 
0.09 0.10 
0.07 0.10 



IMF average of 15 M©, 20 M©, and 25 M© stars: 

Sr 4.402 2.201 0.144 0.072 0.072 

Y 2.844 1.918 0.093 0.063 0.030 

Zr 3.089 1.519 0.101 0.050 0.051 



^Calculated from Eq. 1. 

^f-max - f-min: Spread of fractional production among calculated reaction rates. 



